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Abstract. Using electromagnetic interaction as an example, response transfor- 
mations [L.P. and S.S., Ann. Phys. 323, 1963, 1989 (2008), 324, 600 (2009)] are 
applied to the standard perturbative approach of quantum field theory. This ap- 
proach is rewritten in the form where the place of field propagators is taken by 
the retarded Green function of the field. Unlike in conventional quantum-field- 
theoretical techniques, the concept of space-time propagation of quantized field is 
built into our techniques. 
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1. Introduction 

In this paper we continue our investigation of dynamical response properties of 
quantum systems. In papers [U [2j [3], we introduced response transformations of 
quantum kinematics. In paper |4], response transformations were extended to the key 
technical tool of quantum field theory, Wick's theorem [5j [6l [71 18] . In this paper, we 
rewrite in response representation standard perturbative techniques of quantum field 
theory [71 |9l (TOl jTl] . As a practically important example we consider electromagnetic 
interactions of light and matter. 

Our approach unifies quantum field theory, phase-space techniques and Kubo's 
linear response theory, the latter generalised to nonlinear and stochastic response. 
As was explained in the introduction to a previous paper [1] , this approach is strictly 
subject to the Schwinger-Perel-Keldysh type of techniques jnHTUldTI and does not hold 
(in fact, cannot be even formulated) in the conventional Feynman-Dyson framework 
[7]. We make extensive use of functional methods of quantum field theory; for an 
excellent introduction to these we refer the reader to Vasil'ev's textbook [8|. Functional 
(Hori's [6j) form of Wick's theorem is discussed, e.g., by Vasil'ev and in our paper [4]. 

The paper is structured as follows. In section [51 we take a "bird's eye view" 
at the physical motivation and the key results of this and two forthcomming papers 
[121 [13]. In section |3l we summarise formal definitions. In section |4l we consider 
the broad-band, or nonresonant, case characterised by the real-field-times-real-current 
form of electromagnetic interaction. We apply response transformations to Dyson's 
perturbative techniques. Wick's theorem in its response form (the causal Wick 
theorem [J], for short) emerges naturally in transformed perturbative relations. In the 
broad-band case, these relations stay relatively compact, which allows us to present 
the formal reasoning more or less in detail. The narrow-band, or resonant, case 
considered in section [5] is much bulkier. However, the logic changes little, so we 
only list the key intermediate relations and final results. In section |6l we show that 
the two cases can be seamlessly merged in a single problem. In [Appendix A[ we 
construct a closed diagrammatic solution to the problem of "dressing" of matter by 
electromagnetic interaction. This appendix constitutes a formal "closure" of the paper, 
by associating key dynamical relations with a conventional computational method: 
Wyld-type diagram techniques [2] . In [Appendix Bj we demonstrate that the narrow- 
band case is indeed a resonance approximation to the broad-band case. The former 
may be recovered by dropping the counter-rotating terms in formulae for the latter. 
Although as expected, this result is important for overall consistency. 

2. Bird's eye view of macroscopic quantum electrodynamics 

To allow the reader a better insight into our physical motivation, we start from a 
quick overview of this and two forthcoming papers. A convenient starting point is 
Fermi's famous arrangement [15l [16l [T71 [HI [HI [20] depicted schematically in Fig. 
[1] Two spatially separated atoms are prepared, one in the ground and the other in 
the excited state. The atoms are coupled to quantised electromagnetic field. How 
quickly the state of the second atom would start changing due to the presence of 
the first atom! The italicized reservation is essential, because the state of the second 
atom starts changing immediately due to its interaction with the electromagnetic field. 
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Figure 1. Schematics of Fermi's arrangement. 



Such transients plague any theory with interactions switched on and off — to which 
preparation of atoms in their bare ground and excited states is obviously equivalent. 
One can think of a number of ways how to eliminate this transient. One may look at 
the difference between transients with the first atom present or absent [19] , or improve 
on the concept of state preparation [20 . However, the most consistent — and most 
physical — way is to adhere to the quantum-field-theoretical viewpoint where all 
interactions are switched on adiabatically in the remote past . Formally, this means 
working in the Heisenberg picture, with all physics expressed in terms of Heisenberg 
operators describing field and matter. The price to pay is that state preparation at 
a finite time is not allowed, and we have to look for more realistic models involving 
explicit mechanisms of excitation of atoms. 

Rather than formulating a specific model of Fermi's arrangement, we consider a 
general case of two distinguishable devices coupled to quantised electromagnetic field. 
Such system is described by the generic Hamiltonian, 

N 

H{t)^hJ2uj^aia^ + HdoAt)+Hi{t)- (1) 

By definition, we write all Hamiltonian in the interaction picture. The oscillators are 
represented by the standard bosonic creation and annihilation operators, 

[a^,al,] = SkK', k,k'^1,---,N. (2) 

The nature of the mode index is arbitrary; nothing prevents it from being continuous. 
The mode frequencies are also arbitrary. The interaction Hamiltonian has the standard 
field-times-current form, 

Hi{t) ^ ~ j dxA{x,t)J{x,t). (3) 

The operator of the quantised field t) is an equally standard linear combination 
of the creation and annihilation operators, 

N r-j— 

A{x, t)=Y, J —u^[x)a^e-'^-' + H.c. , (4) 

K— 1 

where Uk{x) are complex mode functions. Variable x comprises all field arguments 
except time. Its meaning, as well as the meaning of the symbol J dx, is problem- 
specific. The Heisenberg (initial) quantum state of the system p factorises into that 
of the oscilators and that of the device. The initial state of all oscillators is vacuum, 

/5= |0)(0| (^Pdcv (5) 




Figure 2. Schematics of interactions in the two-device modeh conventional 
quantum approach. 




Figure 3. Schematics of interactions in the two-device model: field properties 
are expressed in terms of the retarded Green function. 

The device Hamiltonian Hdcv{t), the current operator J{x,t) and the state of the 
device pdev are placeholders for quantum properties of matter. 

In mathematical terms, the Hilbert space of the system is a direct product of 
the field and matter subspaces. The creation and annihilation operators are defined 
in the former, while Hdcvit), J{x,t) and pdov — in the latter. The assumption of 
there being two distinguishable devices is introduced postulating that the matter 
subspace factorises into two subspaces where the triads HdcvA{t), JA{x,t), pdevA 
and Hdevsit), Jsix^t), pdcvB are defined. These operators characterise macroscopic 
components of the device, which we refer to as "device A" and "device B." By 
definition, for operators characterising the composite device we have, 

Hdcv{t) — HdcvA{t) + HdcvB{t)- (6) 
J{x,t) ^ JA{x,t) + JBix,t). (7) 
Pdcv = PdcvA PdcvB- (8) 

Apart from most general assumptions of Hermiticity and positivity all operators here 
are arbitrary. 

Structure of electromagnetic interactions in the outlined model is illustrated 
in Fig. [21 Assumed distinguishability of the devices leads to separation of the 
electromagnetic self-action problems for the devices from the problem of their 
electromagnetic interactior^ For Fermi's arrangement, this means that one may 
separate transients from the interaction of atoms via the field (cf. also endnote [32]). 
Furthermore, both in Fig. [T] and in Fig. [5] the electromagnetic interaction is shown 
as bidirectional. In Fermi's arrangement, single-directedness is achieved by preparing 
the atoms in special quantum states (subject to reservations made in endnote [32 ). 
Can one introduce single-directedness in the general arrangement of Fig. [2P 



X This statement as well as similar statements below are to be formally verified in this and 
forthcoming papers [121 113) 
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Figure 4. Approximations and modifications in the two-device model. Top: 
neglecting macroscopic back-action leads to generalised photodetection theory. 
Bottom: full quantum description of a device is restored by adding c-number 
sources to the model. 

Formally, directedness of electromagnetic interaction is associated with the linear 
response function (also known as retarded Green function) of the free electromagnetic 
field. This quantity emerges if replacing the current operator in ([3]) by a c-number 
source Je{x,t). This decouples matter from interaction. One can then calculate the 
linear response of the field [UEI] (see also [22]). 

A^ix, x',t-t') = ^-^Trrrr^ U = ^^(* - [^(^' ' 

oJe[x, t'j I Jo=0 n 

where A{x, t) is the Heisenberg field operator. The intermediate expression here is 
the definition of Ar, and the last one is Kubo's formula for it [21]. Kubo's proper 
expression contains quantum averaging of the commutator, which we omitted because 
the free-field commutator is a c-number anyway. Kubo's formula may be inverted 
resulting in the wave quantisation relation [23], 

[A{x, t), A{x', t')] = -ih[An{x, x',t-t')~ An{x',x, t' - t)] . (10) 

Commutator (jlOp is all one needs to know about the field in order to construct the 
standard nonstationary perturbation solution to Hamiltonian ([1]). Technical means 
allowing one, firstly, to construct generalised perturbative series without specifying 
the model of matter, and, secondly, to rewrite it (series) in terms of Ar is rather 
involved, but conceptually everything remains fairly simple. The resulting structure 
of electromagnetic interactions is illustrated in Fig. [Sj 

The roles the concept of single-directedness plays in the self-action and in the 
interaction problems are drastically different. In the case of self-action, both ends 
of Ar are attached, so to speak, to the same electron. There is no way to control 
microscopic actions and back-actions associated with emitting and reabsorbing the 
field by matter at a microscopic level. Whereas in case of interaction, the signals 
travelling from device A to device B may be experimentally distinguished from 
those travelling from device B to device A. This assumption concerns not only 
physics but also engineering. Physics implies that the situation is macroscopic (often 
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termed mesoscopic, see endnote [33]): devices are distinguishable and separated by 
macroscopic distances, and all light beams may be controlled. Engineering takes care 
of such problems as, for instance, preventing light reflected off the detector input 
window from affecting the laser source. 

Distinguishability of the incoming and outgoing electromagnetic signals for 
macroscopic electromagnetic devices are put to use in two ways. Dropping the signal 
travelling from device B to device A yields a generalised photodetection theory (Fig.HJ 
top). The dropped incoming wave may be restored by introducing c- number external 
sources into the Hamiltonian (Fig. 21 bottom). It is not immediately clear that c- 
number sources suffice to describe the incoming signal which in arrangement of Fig. [3] 
may be in a quantum state. This statement constitutes a generalisation of Sudarshan's 
optical equivalence theorem |24| to arbitrary interacting systems. 

The said generalisation of Sudarshan's theorem is a particular case of a more 
general statement. Namely, in the formal structure corresponding to Fig. [3J Planck's 
constant turns out to be eliminated from all relations between the quantised field 
and current [23] . Any such relation survives the limit /i — >■ unchanged, and must 
therefore correspond to a particular relation of classical stochastic electrodynamics. 
The most natural expression of such quantum-classical corrrespondences emerges by 
mapping quantum relations into phase-space. 

The order in which the aforesaid formal facts will be established is reversed 
compared to the logic of the above "pictorial argument." The subject of the 
present paper is a solitary device interacting with c-number sources, as illustrated 
in Fig. |4j bottom. Phase-space mappings and quantum-classical correspondences 
will be the subject of paper [12]. Interactions of distinguishable devices, generalised 
photodetection theory and the optical equivalence theorem will be the subject of paper 
[13]. Apart from interaction ([3]), in all three papers we also look at its modification 
under the rotating wave approximation, and at how these two types of interaction 
coexist in a single model. 

3. Generic model of a macroscopic electromagnetic device 

3.1. RWA or no-RWA? 

We consider a structural model of electromagnetic interaction where two sets of 
oscillator modes interact with a "quantum device." The nature of the "device" may 
be arbitrary. For the modes, we assume two different types of coupling: a resonant 
electrical field-time-dipole coupling for one set, and a nonresonant potential-times- 
current coupling for the other. Later this will allow us to treat on an equal footing 
the light (of which the interaction with matter in quantum optics is commonly treated 
under the resonance, or rotating wave, approximation — RWA) and the photocurrent 
(which is a broad-band process and in no way is a subject to the RWA). We cannot 
just assume the RWA because we then lose the low frequency photocurrent modes. 
Nor are we really keen to keep our analyses without the RWA because this would 
essentially disconnect us from the quantum-optical paradigm. (To recognise the 
problem, the reader may try, for example, to reformulate Glauber's photodetection 
theory [25ll26ll22] in terms of real fields rather than analytical signals.) There is no 
simple notational compromise because the number of dynamical variables in the two 
cases differ: one real field versus a pair of conjugates. We therefore see no other choice 
but to develop two versions of the theory in parallel. 
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To refer to physics rather than to formal techniques we shah talk about the broad- 
band (or nonresonant) and the narrow-band (or resonant) cases. Both the response 
transformations and the causal Wick theorems in the two cases differ. The broad-band 
case is a continuation of our papers (and of paper j23j). The narrow-band case 

is an extension of paper [28] , see also [29] . In paper [4] we referred to these cases as to 
the real and nonrelativistic ones, respectively. The change of terminology follows the 
change of motivation: advancing formal techniques in paper [4] versus understanding 
macroscopic quantum electrodynamics in this and forthcoming papers. 



3.2. The Hamiltonian 

Formally, we assume Hamiltonian ([Ij with some amendments. All field oscillators 
are divided in two groups of M and N — M modes, 1 < Af < iV, organised in two 
electromagnetic-field operators. The narrow-band^ or resonant^ field is, 

,.=1 ^^^^ 

E\x,t) = -^Y,^ul{x)aie^^-^-^>)\ 

K—l 

and the broad-band^ or nonresonant, field is defined as. 



N 



A{x,t)^ J^u^{x)a,e~''^-' +II.C.. (12) 

K = M+1 * ^'^ 

In pT]) . Wo is the characteristic optical carrier frequency. Other quantities in (|lll) and 
(|12p have the same meaning as in section [2] For simplicity we assume that both fields 
depend on the same set of arguments. This restriction may be easily lifted. 

For the resonant modes, their frequencies are supposed to occupy a narrow 
band of frequencies centered at loq. That is, E{x, t) and E'^{x, t) are by definition slow 
amplitudes. However, this assumption is only important for physics. Formally, it may 
be disregraded, with the only exception of [Appendix B In the rest of the paper, ljq 
is an arbitrary quantity which may safely be set to zero. 

Electromagnetic interaction in (|T]) is also divided in two, 

Hi{t)=H'l'\t)+H[^\t). (13) 
The narrow-band field interacts with the device according to the resonant Hamiltonian, 

H[^\t) = - j dx[b{x,t)E\x,t) + b{x,t)El{x,t) + D^{x,t)E\x,t)\ H.c. , (14) 

while the broad-band field — according to the nonresonant Hamiltonian, 

H[^\t) = - j dx[j{x,t)A{x,t) + j{x,t)A^{x,t) + Jc{x,t)A{x,t)\. (15) 

The state of the system is given by equation Motivation for introducing c-number 
sources into the interaction was elucidated in section [2] 

Comments in the paragraph preceding equation ([6]) hold, except there is an 
additional operator D{x,t) characterising matter. Factorisation of the matter 
subspace will not be assumed till paper [13], so that equations ([6]), (O and ^ do 
not apply. 
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Fields, currents and dipoles in equations (1151) and (I12[) are interaction-picture 
operators. Their Heisenberg counterparts will be denoted by calligraphic letters as 
£{x,t), A{x,t), T>{x,t), and J{x,t). 

Hamiltonian ([T|) is a placeholder for all conceivable cases of electromagnetic 
interaction, from a single mode (with M = 1) to relativistic quantum fields (with 
M = and N — >■ oo). From our perspective, this Hamiltonian is a structural model 
of a quantum-optical experiment involving photodetection. Operators £{x,t) and 
T>{x,t) describe optical interactions, while A{x,t) and j{x,t) express photovoltage 
and photocurrent. Without specifications, Hamiltonian ([1]) is a structural model of 
photodetector. By breaking the device into distinguishable components, we recover a 
structural model of the source-detector interaction, etc. 

3.3. Operator orderings 

Here we summarise definitions of operator orderings used in the paper. All definitions 
imply bosonic operators. Of importance to us are the normal and time orderings. The 
normal ordering puts a^'s occuring in equation ([1]) to the left of a's, 

■anO^^r. ^ ■a^^,a>,-. ^ a\,ai^, (16) 

etc. It is extended to the field operators ((TT|) . by linearity. The time-orderings 
T+ and T_ put operators in the order of decreasing and increasing time arguments, 
respectively, 

T+Xi{ti)X2{t2) ■ ■ ■ Xm{tm) 

= 0itl-t2)0{t2-t3)---0{t„,^i-t„,)Xi{tl)X2{t2)---X,nit^), (17) 

perms 

T-Xi{tx)X2{t2)---Xn{t^) 

= H 0{tl-t2)0{t2-t^)---e{tra-l-tra)X„,{tra)---X2{t2)Xi{ti). (18) 

perms 

The sums are over all permutations of times ti^ • • • ^tm- Of use is the property, 

[r±A'i(ti)i'2(t2)---i'™(t™)]^ = r^^/(ti)^i(i2)---i'i(im)- (19) 

It expresses the simple fact that Hermitian conjugation inverts the order of factors 
turning a r+-ordered product into a T_-ordered one, and vice versa. We note in 
passing that definitions like P7)) . (IT51) may cause mathematical problems, see the 
concluding remark in section IIIB of paper [4 . 

In formulae, the Tj- -orderings mostly occur as double time ordered structures 
T_ ■ • • r_|_ • • •. Rather than visually keeping the operators under the T±-orderings, one 
marks the operators with the ± indices and allows them to commute freely: 

T_i'i(ti) • • • x^{t^) T+yi{t[) ■ ■ ■ 3>„(4) 

= TcXi^iti) ■ ■■Xrn-{tm)yi+{t[) ■ ' •i'n+CC) (20) 
= Tcyi+{t[) ■ ■ ■ y„+{t'JXi-{tl) ■ ■ ■ Xrn-{tm), 

etc. The ± indices serve only for ordering purposes and otherwise should be 
disregarded. Using Tc results in a reduction of formulae in bulk which may be truly 
dramatic (compare, e.g., expressions for response functions in [2] and in [3]). In this 
paper we mostly use the Tc-ordering and revert to the double time ordering only 
where necessary. 
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3.4- Quantum Green functions 

Quantized fields E{x,t) and Aix^t) enter all formulae through their Green functions. 
In either case, we define the retarded Green function and two contractions: the 
Feynman propagator, and one more kernel which emerges as a propagator in Perel- 
Keldysh's diagram techniques (and which, strictly speaking, is not a Green function). 
For the resonant field, these quantities read, respectively, 

Gk{x,x' ,t-t') = lea - t')(Q\\E(x,t), EUx' ,t')]\Q), (21) 

h ' ' 

GAx,x\t-t') = l(Q\T+Eix,t)EUx' ,t')\Q), (22) 
ft ^ ' ' ' 

G'^+Xx.x'^t-t') ^ ^-{0\E{x,t)E\x',t')\0). (23) 
For the nonresonant field, they are, 

AK{x,x',t-t') = ^-e{t-t'){Q\[A{x,t),A{x',t')\\Q), (24) 

AF{x,x\t-t') = ^{Q\T+A{x,t)A{x' ,t')\Q), (25) 

A'^+\x,x\t-t') = ^{0\A{x,t)A{x',t')\0). (26) 

in (|22p and ([SS]) is the standard time-ordering defined by equation (|17|). Using 
(jlip . (IT^ we find the explicit formulae, 

G'-^^x.x ,t-t) =«2^ ^ , 



K = l 

N 



(27) 



2wk 

K = M + 1 



Other kernels follow from the relations, 

Gn{x, x', t-t') = Gf{x, x', t^t')^ e{t - t')G'^+^ [x, a;', t - t'), 

AR(a;, x\ t-t')^ e{t - t') {^x, x\ t-t')- A(+)(a;', x, t' - t)] , (28) 

Af(x, x',t~t')^ e{t - t')A'+^ (x, x\ t-t')+ e{t' - t)A(+) (x', X, t' - i). 

Definitions of the retarded Green functions (PT|) and ([Ml) are Kubo's formulae for 
linear response functions [H]. Definition of Ar was discussed in section[51 cf. equation 
([9]) and comments thereon. Gr is defined in a similar manner, setting D{x,t) = in 
Hamiltonian and then defining the linear response of the resonant field by the 
formula. 

Equation (PT|) is Kubo's formula for Gr. For more details on the linear response theory 
for free bosonic fields see papers [Tl[22]. Commutators in ((2T|) and (j24p are c- numbers, 
so that the averaging may be dropped (in other words, response of a linear system 
does not depend on its state). In the other four definitions, this averaging matters. 
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3.5. Response transformation of contractions 

All our results hinge on one observation: the propagators (contractions) may be 
expressed by the retarded Green functions [U [H [531 HH]- In the narrow-band case, 
this connection readily follows from (l28l) . 

GF(x,x',t-t') ^Gn{x,x't-t'), 

G'^+\x,x',t-t')^GB.{t-t')-G*^{x,x',t' ^t). 

In the broad-band case, such connection is more involved. It happens to be an integral 
transformation [J HI [5S] , 

Af{x, x',t-t')^ A^+\x, x', t-t') + a\^\x\ X, t'~t), ^^^^ 

A(+) {x, x\t-t')= a|,+' {x, x',t-t')- A^^'> {x', X, t' ~t). 

The symbols denote separation of the frequency-positive and negative parts of 
functions, 

fit) = f(+\t) + f(-\t), 



J-oo 



(32) 

= / dte'^^'fit). 



H-oo 

iujt , 



The frequency-positive and negative parts are always defined with respect to 
native arguments of functions; so, A'^\x' , x , t' — t) in ([5T|) are deciphered as 
A^^(x',x,i)|j^^,_j. For more details on this operation see paper P], appendix A. 

Note that notation for A(+) and G^^^ is consistent with them being purely frequency- 
positive. 

Equations pip may be traced down to the wave quantisation formula (|10p , while 
(|30p — to a similar formula for the narrow-band field, 

[E{x,t),E''{x',t')] ^ -ih[Gn{x,x',t-t') -G'^{x',x,t' -t)]. (33) 

Note that Ar is real, while Gr is complex. For more detals see papers pi HI 128). 

3.6. Condensed notation 

To keep the bulk of formulae under the lid and make their structure more transparent, 
we make extensive use of condensed notation, 



fg^ J dxdtf{x,t)g{x,t), (34) 

fKg^ J dxdx'dtdt'f{x,t)K{x,x',t-t')g{x',t'), (35) 

{Kg){x,t)^ J dx'dt'K{x,x',t-t')g{x',t'), (36) 

{fK){x, t)= dx'dt'g{x', t')K{x', X, t' - t), (37) 



where f{x,t) and g{x,t) are c-number or q-number functions, and K{x,x' ,t — t') is a 
c-number kernel. The "products" fg and fKg denote scalars, while Kg and fK — 
functions (fields). 
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4. Oscillators interacting with quantised current (the broad-band case) 

4.1. The closed-time-loop formalism 

Perturbative formulae for Hamiltonian ([T]) grow insufferably bulky. We therefore start 
from the relatively compact broad-band case. Formally, this corresponds to setting 
M = in equations ©-([121). 

To maintain connection with papers [TJ [2J [31 [221 [2E] > we solve for the characteristic 
functional of products of the Heisenberg field and current operators, A{x,t) and 
J'{x,t), ordered in the Schwinger-Perel-Keldysh closed-time-loop style [TUl [TTj . 

2(f7+,?7-,C+,C-l-4o, Jc) 

= (T_ exp ( - - iC-j) T+ exp (ir^+A + iC,+j)) (38) 

= {Tc exp {iri+A+ - irj^A- + iC,+J+ - iC-J-)), 

where ri±{x, t) and C±{^j t) ^-re four independent complex functions. We use condensed 
notation p4p . The averaging in ([55)1 is over the initial (Heisenberg) state of the system, 

(•••)=Trp(...), (39) 

where p is given by equation ^ and the ellipsis stands for an arbitrary operator. 
Definitions of operator orderings were summarised in section 13.31 

In this paper, we denote various characteristic functionals of Tc-ordered operator 
averages as S, while the same functionals in causal variables will be denoted as 3>. 
This agrees with notational conventions of paper [21] . Notation used in papers [1] [31 [3] 
is recovered replacing S ^ $ and $ ^r. The reason for the change of notation 
compared to [Tl[31|3] is obvious: we attach quite a number of indices to the functionals, 
and an additional index would be a nuisance. 



4-2. Dyson's perturbative technique 

As in papers [23l [28] we employ Dyson's standard perturbative techniques [7]. We 
assume that the reader is familiar with the concept of S-matrix and its representation 
as a T-exponent, 



S = T+ exp 



dtH{t) 



(40) 



where H{t) is the interaction Hamiltonian in the interaction picture. As in papers 
[nisi [31 [5, omitted integration limits imply the maximal possible area of integration: 
the whole time axis, the whole space, etc. 

The key relation for us is the formula expressing time-ordered products of 
Heisenberg operators, Xi{ti) , ■ ■ ■ , Xm{tm) , as time-ordered products of the same 
operators in the interaction picture, Xi{ti), • • • , Xm{tm), 

T+Xi{ti) ■ ■ ■ X^{tra) = S^T+SXi{ti) ■ ■ ■ X^{tm). (41) 

Detailed derivation of this formula may be found, for instance, in Schweber's textbook 
[7]. Uncharacteristically for ever-so-thorough Schweber, derivation in ,7, disregards the 
iS^ factor; for the necessary amendments see our paper . section 5.2. Taking the 
adjoint of equation (|^Tt and using equation we also obtain, 

T^Xi{t^) ■ ■ ■ X^{t^) ^ [T^S^X^{ti)---X^{t^)\S. (42) 

Square brackets mark the region to which the ordering applies. 
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Equations (pUjl - dl^ imply that the S-matrix is regarded a functional of the 
operators A{x,t) and J(x,t). The reader should keep this in mind to avoid confusion 
and misunderstanding. For details see Vasil'ev [S]. 

In the case at hand, the S-matrix and its inverse read, 



S = T-f exp 
5t ^ T_ exp 



■{AJ + AcJ + AJc 



Uaj + AcJ + AJc 



(43) 
(44) 



where A{x, t) and J{x, t) are the field and current operators in the interaction picture. 
The former is given by (jl2[) . while the latter is just assumed to be known. We continue 
using condensed notation (l34l) . Applying equations (|4T|) . p2|) we obtain, 

T+ exp {ir]+A + iC+J) 

= 5^r+ exp {irjj^A + «C+ J) exp 
T_exp(-i?7_yl-<_j) 

T_ exp ( — irj^A — i(^-J) exp — — 



-(AJ + AcJ + AJc 

h 



{Aj + A^J + AJc) 



(45) 



(46) 



Unlike in equations (|4T|) . (|42|) . in equations p3)) -(|46 | the distinction between operators 
<S and iS^ and their representation as time-ordered operator functionals is explicitly 
maintained. When substituting equations (j45|). (j46|) into (p8| . the and 5 factors 
cancel each other, 

ss^ = i. 



(47) 



and we have, 

'^(v+,v-,C+,C-\Jc,Ac 



■{Aj + A^J + AJc) 



■{Aj + AeJ + AJc) 



T_ exp ( — irj^A — iC-J) exp 

xT+ exp {irj-f-A + iC+J) exp 
The same in terms of Tp-ordcring reads, 

'^{v+,V-X+X-\Jc,Ac) = / Tc exp (i?7+i+ - i?7_i_ -iC-J- 



(48) 



X exp 



■ {A+j+ + aJ+ + A+j^) 



exp 



(49) 



Grouping of arguments in this formula is an allusion to its "ancestors:" equations 
(EHl), (SI and (HI. 



4-3. Wick's theorem and elimination of field operators 

In the standard techniques of quantum field theory [71 IB] , the next step is applying 
Wick's theorem [5] to equation (P^)) . Wick's theorem allows one to rewrite time- 
ordered operator structures in a normally-ordered form (for definitions of operator 
orderings see section 13. 3p . This comes especially handy if the initial state of the 
system is vacuum. 
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The interaction-picture operators t) and J{x, t) are defined in two orthogonal 
subspaces and may be manipulated independently. Due to assumed factorisation of 
the density matrix, this equally applies to quantum averaging. We therefore apply 
Wick's theorem only to the field operators in (HHl) . keeping the current operators "as 
are." 

We employ the functional (Hori's [6]) form of Wick's theorem (the Hori-Wick 
theorem, for short). It is discussed by Vasil'ev [8] and in our paper [4]. The directly 
applicable case is that of real field, see paper [4] , section VIII. The Hori-Wick theorem 
for a real field reads, 

Tc^(i+,i_) = :J-A,(i,i):, 



Fn{A+,A^) = expZc 



(50) 



^5A+' 5A. 

where A±{x, t) are auxiliary c-number functional arguments, and J-'(-, •) is an arbitrary 
c-number functional. The reordering form Zc{-^ ■) is a functional bilinear form where 
the kernels are contractions, (with f±{x,t) being auxiliary functional arguments) 



Zc{f+J-) = -f /+Af/+ + f /-a;/_ - ^fi/-A(+)/+. 

We use notation ([35]). Differentiations are carried out using the relations, 
SA+{x,t) SA^{x,t) 



6A+{x',t') 6A_{x',t') 
SA+{x,t) 5A^{x,t) 



S{x-x')S{t-t'), 



(51) 



(52) 



6A-{x',t') SA+{x',t') 

Note that under the normal ordering in equation ([50)) the ± indices of the operators 
are dropped as irrelevant. 
For any functional -Fjv('i 

(0|:J-w(i,i):|0) = J-Ar(0,0), 



(53) 



and we find, 

'^{v+,V-,C+,C-\Jc,Ac) =expZc 
\. J, 

X exp 



f '5 



+ - A+-^hJ_ + - A 



\SA+' 6A^ 



xTrpdcvrc-exp 



Z J4- 



A+ + A^ 



+ C-+ 



iJ- 



A-+A^ 



(54) 



A±=0 



The necessary and sufficient information about the device is thus collected in the 
functional, 



"dcv(C+, C-) = TrpdcvTc exp (iC+J+ - <_ J_). 



(55) 



With it we have. 



^{v+,V-,C+,C-\Jc,Ac) =expZc 



X exp 



A+-i[v 



/_5 S_ 

[SA^'SA 

Jc 

h 



x = 



'dcv 



A, 



Ao , A_ 
— + (+,— 



A, 



(56) 
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Response transformations revisited 

We are interested not in functional S as such but in its response transformation [TJIHIS], 
(cf. the closing remark in section HTTj) 

^('y^CljcaclJc, Ac) = S(77+,77^,C+,C-|-'c, Ae)|^.^., (57) 

where c.v. (short for causal variables [Jl [2H]) refers to the response substitutions 

[3 man], 

ic{x,t) 



'n±{x,t) 



C±{x,t) = ^^±C^^\x,t). 



(58) 
(59) 



The frequency-positive and negative parts are defined by equation ([32]). Formulae for 
the causal variables are, 

r]{x,t) ^ri+{x,t) -r]^{x,t), j,{x,t) ^ n[r/:+\x,t) + T^^r\x,t)], (60) 

ax,t)^C+{x,t)~C^{x,t), a,{x,t) = h[C^+\x,t)+C^-\x,t)]. (61) 

Equations (|58|) . ([SQ)) and (|60|) . (|6ip are mutually inverse, showing that these formulae 
constitute a genuine change of functional variables. Response substitutions are 
inherently related to equations ([5T|). For details see our paper pj. 



4-5. The causal Wick theorem revisited 

For the two exponents in (|56p . the laws of transformation were worked out in paper [4], 
section VIIIB. With = 1 and Jc{x,t) = 0, equation (|56|) reduces to the "test-case 
formula" (104) in paper [3]. To rewrite that relation in causal variables we applied 
additional response substitution, 

A±{x,t) = a'^{x,t)±hC^'^'>{x,t). (62) 

Variables a'^{x,t), C,'{x,t) are given by the formulae, 

A+{x,t) - A^{x,t) 



C{x,t) = 



<(a;, t) = A^^^ (x, t) + AL"^ {x, t) 



(63) 



Equations are a rescaled version of ([5^ . (pTj) : they turn into the latter if 

we define, 

A±{x,t) 



C'±[x,t) = 



Primes keep variables (|63p distinct from (|6T|) . 
In causal variables, 



I ?7+ + 



A+ - M '7- 



ft 



and 



/ 5 



\5A+' 5A^ 
where we use condensed notation ([34|) . (|35|) . 



(64) 



= ^ac + C'(jc + Jc), (65) 



(66) 
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4.6. Closed perturbative formula in causal variables 
As to the functional Sjj^^, we define, 

<&dcv(C|ac) =Si„,(C+,C-)|,.v., (67) 
where c.v. refers to substitution (pT|) . Combining ([5T|) and (p^ it is easy to show that, 

^M^i^l±^l±i^ + C±(x,0 

= M^,t) + a^{x,t) + a'^{x,t) ^ ^(^) ^ ^,(^) 

hence 



77I 



+ C+ /°^^' +C-) -$Lv(C + C'|Ae + ac + aO. 

/ r. V. 



(69) 



Putting equations ([65)) . ([66j) and (j69|) together we obtain the relation sought, 

/ 6 6 
*(?7,C|jo,ac| Jc^e) = exp f - «^^R^ 

X exp [iTja'^ + iC'{jc + Jc)]$dcv(C + C'l^c + (70) 
The generahsed consistency condition [3J [3] , 

$(?7, C|jo, ac| Jo, ^c) = ^{V, (\jc + Jc, Oc + Ao|0, 0) EE C|io + Jc, Ac + A), (71) 

is evident in (ITni) . 

In what follows, we set ^e(a;,i) and Je{x,t) to zero. Relations with sources, 
which are important for physics, may always be restored replacing ae{x, t) ae{x, t) + 
Ac{x,t) and jc{x,t) jc{x,t) + Jc{x,t). 



4-7. Dressing the current 

Unlike in the linear case [4], finding a closed solution to ([70| is certainly impossible. 
One can, however, reduce (170]) to the functional, 

$dcv(C|ac) = $(0,C|0,ae) = $(0,C|0,ae|0,0), (72) 

cf. equation ([7T|) . This functional expresses the properties of the Heisenberg current 
operator (physical, or "dressed," current). A perturbative formula for $dcv is given 
by equation ([70| with ri{x,t) = jc{x,t) ~ Ao(a;,t) = Jc{x,t) = 0, 

$dcv(C|ac) - exp (^-z^AR,^^$Ji,,(C|ac). (73) 

Functional expresses the properties of the interaction-pictire current operator 

("bare" current). Equation ([75]) thus expresses "dressing" of the current operator by 
the electromagnetic interaction. 

Finding the dressed current is the hard part of the problem. Vast majority 
of technical developments in quantum field theory, condensed matter physics and 
quantum optics are in essense attempts to approximate equation ([73|) or its 
counterparts. Connections between equation ()73p and diagram techniques are 
elucidated by Vasil'ev (8j. In [Appendix A| we construct a formal solution to ([73]) 
in terms of a diagram expansion. That (I75|) may be used as an alternative entry point 
to phase-space techniques was shown in paper |28) . 



Quantum electrodynamics in response representation 
4-8. Solution in terms of dressed current 



16 



Assuming that equation ([75)1 is solved, full electromagnetic properties of the device 
follow with ease. We employ the formula, 



exp I f-^ j exp [igh) ~ exp [igh) exp 



S_ 



ih 



(74) 



where f{x,t), g{x,t) and h{x,t) are c-number functions. It can be verified, e.g., 
expanding all exponents in Taylor series. With it we can pull exp [irja'^ + iCjc) from 
under the differentiation resulting in, 

C|jo, flo) = exp {iija'^ + iCjc) 



X exp 



Setting C'(x,t) and a'^{x,t) to zero we find. 



exp 



5 

(5ao 



*dcv(C|ac). 



(75) 



(76) 



Bilinearity of the form in the exponent allows us to write. 



exp 



SOo. 



irj Ar 



5_ 



= exp (ii^Aujc) exp ( — ArJc ) exp ( V^Rj^ ) exp ( - i—^^J? 



'6C 



(77) 



The last factor here is the dressing operator, cf. equation ((73|). The second and the 
third ones are functional shift operators. Applying the dressing and the two shifts to 



we obtain. 



^{riX\jc,ac\Jc, Ac) 

= exp [i77AR(jc + Je)]*dcv(C + ?7^R|ac + + Au{jc + Jc)), 



(78) 



where we use abbreviated notation of equations and ([57| . In equation ([75)1 the 
external sources are restored. Without the device (^dov = 1) and with Jc{x,t) = we 
recover the test-case formula of paper [3] , 

*vac {v\3c) = exp (irjAjxjc) . (79) 
The physical content of equations (1751) and ([75]) is a subject of forthcoming papers. 



5. Oscillators interacting v^rith quantised dipole (the narrow-band case) 
5.1. Closed-time-loop formalism and perturbation theory 

We now consider the second generic case: a quantised dipole interacting with a set 
of oscillators under the RWA. Formally, one sets M = in equations dTTI) . (fT2)) . 
To a large extent, amendments to the broad-band case reduce to considering two 
fields, A{x,t) £{x,t),£^x,t), and two "currents," J{x,t) 'D{x,t),'D\x,t), with 
£{x,t) coupled to &{x,t) and £^{x,t) to 'D{x,t). We remind that £{x,t), V{x,t) are 
Heisenberg counterparts of the interaction-picture operators E{x,t), I){x,t); E{x,t) 
is given by equation (fTTj) and D{x,t) is "just assumed known." Consequently all 
arguments in characteristic functionals also double, ri±{x,t) — /i±(a;, i), /I±(a;, i), 
C,±{x,t) i^±{x,t),D±{x,t), etc. The basic object of the theory, the characteristic 
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functional of the Tc-ordered products of the Heisenberg field and dipole operators, 
reads, (cf. the closing remark in section H?lT) 

Tc exp (ip,^£+ + i^j^£\ — iji^E- — ifi-£l_ 



+iD+V+ + iiy+vl - iP^V^ - iv^'b^_)y (80) 
We again employ condensed notation p4p . The S-matrix and its adjoint now are, 

S = T+ exp 1^ [£:L>t + E^b^ + EDI + D + Elt) + E'^ D^] |, (81) 

5t = r_ exp I - [ED^ + E^b^ + EDI + E'' D + E*b + E'^b^]^, 
By the same means as equation (|49p was obtained we find a perturbative formula, 

'^{^i■+,^^+,f|-,^^-,v+,v+,v-,v-\Dc,Dl,Ec,El) 

Tc exp (ifi^E^ + i/i+i?^ — ijl-E- — iii-El_ 

+iD+b+ + iv+bl - iD^b^ - iv^b^_) 



(82) 



X exp <{ ^ [E+bl + E,bl + E+bi + Elb+ + E;b+ + E^b^ 



X exp I - ^ [E^bl + E^bl + E^bl + # £>_ + E*b^ + Elbc]^j. (83) 

The applicable version of the Hori-Wick theorem is that for the complex field in the 
nonrelativistic case (cf. paper 0], eq. (23)): 

TcT{E+,eI,E^,eI) = ■.Fn{E,E\E,E'^):, 

{ 5 5 5 5 \ - - (84) 

p.(^;„ i^,, ^;_) = exp (^_, ^, ^ j^(i^„ i^,, i^_), 

where E±{x, t), E±{x, t) are four auxiliary functional arguments. The reordering form 
Zc reads, (in notation (pSt ) 

Zc{f+, f-J-J-) = -thf+GFf+ + thf^G*pf- - ihf^G^+^f+. (85) 
f±{x,t), f±{x,t) are auxiliary functional arguments. Using that 

{0\:Tn{E,E\E,E''):\0) = J-Ar(0, 0, 0, 0), (86) 
we find the narrow-band counterpart of equation (j56p . 

/- - - - I * *N / 6 S S 6 

'^{tJ-+,l^+,tJ-,y'-,v+,v+,v-,v^\Dc,D^,Ec,E^) = expZ, 

X exp 



5E+ ' 5E+ ' 5E^ ' 5E. 



^[^^ + ^jE++^(^^+ + ^jE+-^(^^l^ + ^jE^-^(^^i^ + ^]E^ 
E*+E+ , E^ + E+ _ , E*+E_ , E^ + E_ 



E±=E±=0 

(87) 

where 

■^dcv^P^^ = TrpdovTc exp {i'^+b+ + iv^b\^ — iD^b^ — iv^b'^_ 
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5.2. Plain versus duplicate phase space 

In paper j^S], techniques similar to those developed here were used to derive various 
generalisations of the positive-P representation of quantum optics (see [53] and 
references therein). The positive-P as well as its generalisations assume doubling 
of the classical phase space. This doubling is also present in response transformations 
for the nonresonant case, derived in paper [4], section IV. They were deduced for the 
"test-case formula" which is nothing but equation (|87l) with S'^"^ = 1 and De{x, t) = 
(i.e., for a free field in a vacuum state). Direct calculation (see paper [J, section VC) 
gives, 

Svac(/i+,M+,Ai-,A*-) = cxpZc(i/I+,iA*+, -«A*-)- (89) 
Response substitutions in the "test-case formula" read, 

^.+ = -, f^+^P + j, -^--j, ^'-=^' + J^ (9o) 

E+^e'^, E+^hv' + e^, E^^e'^, E^ = hv' + e'^, 

where the arguments were dropped for brevity. Referring to these substitutions (or 
their suitable subsets) as "c.v." we find, 

Zc{in+,in+, -ifj.-, ~iJi-)\ ^ = ip-Gndc - inC^d^, (91) 



cf. equation ([35|1. Furthermore, 
5 5 5 5 



Zc 
and. 



5E+ ' 5E+ ' 5E_ ' 5E. 



E{x,t) = J dx'dt'GR{x,x' ,t - t')dc{x' ,t'), 
E{x,t)= I dx'dt'G^{x,x\t-t')d^{x',t'), 



exp {ip.+E+ + in+E+ - iJi^E^ - i^i^EJ) 

= exp {ijle^ + iv dc — i/ie^, — ii^'dc) . (93) 

Doubling of the phase space is best seen in equation ((9T|) which expresses two 
independent fields. 



(94) 



emitted by two independent sources, do{x,t) and dc{x,t). Such formal structure 
comes very handy if the goal of analyses is indeed derivation of generalised positive- 
P representations as in |28) . This and forthcoming papers utilize a version of P- 
representation |30j rather than positive-P, so that doubling of auxiliary variables is 
better be avoided. 

To suppress doubling of the phase space we impose conditions on the causal 
variables, 

p,{x,t) ^ fj.*(x,t), dc{x,t) ^ dl{x,t), 

D'{x,t)^v'*{x,t), e'^{x,t) = e':{x,t). ^ ' 

To preserve response substitutions (^0)) we also have to impose conditions on the 
auxiliary variables, 

ii±{x,t) = ^i*^{x,t), E±{x,t) = E*^{x,t). (96) 
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Conditions (j95p . (j96p are by definition consistent with response substitutions and do 
not interfere with above analyses, nor with analyses in paper [3]. Wick's theorem ((M)) 
equally holds if independent functional derivatives by E±{x,t), E±{x,t) are replaced 
by pairs of complex-conjugate derivatives. All derivations in the above and in paper 
[4] remain applicable. 

Noteworthy is that effective doubling of the phase space is also present in the 
broad-band case. Indeed, if we keep ri±{x,t) real, je{x,t) comes out complex, whereas 
a classical source je{x, t) in ([79|) should be real. To keep the auxiliary variable jc{x, t) 
real, we could introduce the condition, 

r;_(x,t) =77;(x,t). (97) 

However, this makes ri{x, t) purely imaginary, so that analytical extension to real 
ri{x, t) is needed. We leave the question of mathematical rigour behind response 
transformations open for discussion. 



5.3. Response transformation and solution in terms oj dressed dipoles 





dc 
^ ft' 




* 1 

= /i 4 


d* 

' ft' 






Be 
^ ft' 






e* 

ft ' 


e* Co 

h ft 


E+ 


= <, 


E+ 


= hv'* 


+ <, 


E-=e';, E^^hu' + e'^ 



Based on the above, we postulate response substitutions in the narrow-band case to 
be, (again dropping the arguments) 



(98) 



The first and the third lines here are substitutions ([50]) under conditions (^5]). ([M)) . 
The second line is modelled on the first one; it imposes one more condition on auxiliary 
variables, 

v±{x,t) = iy^{x,t). (99) 

Similarly to the broad-band case we define the functionals, (cf. the closing remark in 
section 14.11) 

'S>{^l,^l*,v, v*\dc,dl,ec,el\Dc,Dl,Ec,E*) 

= s(/i+,^+,/2_,^_, iy_|£;c,-B*,i:'c,£'o) lev.' 

$dcv(i^, y*\ec, el) = $(0, 0, V, 0, e^, e:|0, 0, 0, 0), 

with "c.v." referring to (suitable subsets of) equations ((98|) . Rewriting (|87|) in causal 
variables we have, 



(100) 



^[^i,^.* ,v,v*\dc,dl,ec,el\Dc,Dl,Ec,E*) = eyiY> ^ .^^,„„,^^,_^ • ^ 
X exp [i^J*e'^ - i^JLe'* + iv'*{Dc + 4) - w' [Dl + d*)] 

x<^\c,{u + v\ V* + v'*\Ec ee + e^, + e*, + e^*) . (101) 

Confining this to dipoles we find the dressing relation. 



dcv Ky-i ^ ) 



exp 



5e':''5v' 



(102) 
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Using equation ([71| and proceeding in close similarity to equations (|7S 1) " ([75| we find 
the solution in terms of dressed dipoles, 

a**, ly, i^*|dc, die,, el \D,, D:,E,, E;) = exp [iijl*G^{D, + 4) - *AiGJ^(C: + <)] 
x$dcv(j^ + mGr, j^* + m*Gr|£;o + Go + GR(cfo + -Do), + < + GR(cf;; +DI)). 

(103) 

The physical content of equations (|102p and (I103P is a subject of forthcoming papers. 

6. Putting the broad-band and the narrow-band case s together 

There is no difficulty whatsoever in unifying manipulations as in sections 3] and [5] under 
the common umbrella of the model of section [5?^ for arbitrary M and N, 1 < M < N. 
These manipulations apply to different sets of modes, and do not interfere if united 
in a single problem. Formulae for the generic model of section 13.21 may therefore be 
written straightway, by merging the corresponding formulae for the broad-band and 
narrow-band case s. So, the characteristic functional of the Tc-ordered products of the 
Heisenberg field, photocurrent and dipole operators, A{x, t),£{x, t), J{x, t),'D{x, t), is 
a merger of equations (pS)) and ([50)) . (cf. the closing remark in section HTTj) 

'^{v, f^, M*, C, V* I jc, 4, Oo, Cc, e* I Jc, De, Dl,Ac, E^, E*) 

Tc exp (iri+A+ - irj^A- + iC,+J+ - K-J-) 



X exp + i^j^£\_ — i^^S^ — i^^Sl^ 



Mv+V++iv+V\-iv^V^-iv^'D^_))\ , (104) 

/ I c.v. 

where c.v. refers to merger of the broad-band and narrow-band response substitutions 
([55]) . ([55]) and (IM)) . The properties of dressed photocurrents and dipoles are expressed 
by the functional, 

$dev(C, V, V* |ae, ee, = $(0, 0, 0, C, V: v*\Q, 0, 0, Co, e,, el\0, 0, 0, 0, 0, 0). (105) 
The formula reducing ([1041) to ([lOSp is constructed merging equations ([78| and ([103p . 
^, /X*, C, V, V* I jo, 4, 4; «c, Go, e* |4, -De, i ^c, -Bo, -Be ) 
= exp [i?7AR(je + 4) + i/U*GR(4 + -Do) - i^iG*^{dl + D^)] 

V +/i*GR|^ext tot,-E'oxt tot,-E'ext tot), (106) 

where 

^oxt tot (2:, t) = ac(x, t) + Ac(x, t) 



(107) 



(ia;'di'AR(a;, x' ,t — t') [jo (2^', t') -f J^x' , t')] , 
-Eoxt tot(2:,i) = ec(a:;,t) + Ecix,t) 

+ J dx'dt'Gn{x,x' ,t ~ t')[d^{x' ,t') + D^{x' ,t')]. 
Finally, the dressing relation is a merger of ([73p and (|102l) 
$dov(C, |«c, ee, <) = exp [ - ^^Ar A _ ^^GR^ z ^g^^ 

x*dov(C,'^,'^*|ac,eo,e*), (108) 
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with the formula for unifying ([35]), dHT]), dSS]), and piTOl) . 

'J'dev(C,i^:'^*|ac,eo,e*) ^ Trpdcv Jc exp (iC+J+ - iC-^-) 

xexp{iiy+f)+ + iiy+D\_~iD_D_ -iiy_Dl)\ . (109) 



I c.v. 



7. Conclusion and outlook 

General structure of electromagnetic interactions in response picture was elucidated. 
The astonishing feature of perturbative relations thus obtained is that they lack 
Planck's constant. It is hidden in functionals $^1^^ and $dcv and should reappear were 
equations of motion for matter considered explicitly. However, for electromagnetic 
interactions of "dressed" devices Planck's constant is irrelevant. Implications of this 
point will be subject of forthcoming papers [TH [T3] ■ 
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Appendix A. Diagram techniques 

This appendix constitutes a formal closure of the paper, by associating dressing 
relations with a common computational method: Wyld-style diagram techniques |14j . 
It assumes a reader versed in functional techniques of quantum field theory [8]. We 
consider the broad-band case. Our goal is to construct a generalised diagram expansion 
for the dressing relation ([73|) . To this end, the functional is represented in terms 
of bare response cumulants Q'-"^'"\ 



*dov(C|ac) = exp 



E— — / dti--- dUidt[ •••dt^^ 
m\n\ I 



m.n—0 



xQ(™'")(ii, ■■■,t^\t[,--- C)C(ti) • • • C{tMt[) ■ • • ae(C) 



(A.1) 



For simplicity we drop all field arguments except time. The quantities Q^™'") express 
properties of the bare current. For example, in spinor electrodynamics, they are found 
applying response transformations |H [3] to fermionic loops (implying that (lA.ll) is 
properly generalised to the 4- vector density of current). Similarly to functional 
4>dov is represented in terms of dressed response cumulants 



*dov(C|ac) = exp 



— ;— r dti--- dtmdt[ ■■■dt'^ 
'-^ m\n\ J 



m,n— 



X Q("'") (ii, ■ • • , • • • C)C(ii) ■ • • C(im)ae(t'i) • • • ac(C) 



(A.2) 



The quantities Q(™'") express properties of the dressed (observable, physical) current. 
Substituting (jA.ip into (|73l) . expanding all exponents in power series and comparing 
the result to (|A.2p we recover the diagram series for the dressed cumulants. 
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For the sake of argument, we assume that only three cumulants, Q^^'^\ Q^^'^) 
and Q(2,2)^ 

are nonzero. In the diagram expansion, they serve as generaUsed vertices. 



Q(^''Ht\t') = {t'^^ /}, 



Q^'''\t\t[,t'„t's) = { K 




(A.3) 



Q^'''\tut2\t[A) 



t'2 



X ■■ 



Vertices are distinguished by the number of incoming and outgoing lines; their 
orientation and other details arc of no importance;. Curly brackets isolate graphical 
elements visually. Wc also introduce a graphical notation for the causal propagator, 

AR(t - t') = {f * t} . (A.4) 

Then, for example. 



(A.5) 




Q(^'^\t\t[A,t',) = I 



Q^^'^\h,t2\t[,t',) ^ < 



\ I 



> + 



(A.6) 



(A.7) 




t]^ j dhdt2Q^^^^Ht\h)An{ti - t2)Q'^'^'\h\t'), (A.8 



t } — 2 I dtidt2dtzdt4,dt5dtQ 



xQ^^'^\t,ti\t2,t3)AR{t2-t4)AR{t3-t5)AR{te-ti)Q^^'^\t4,t5\t',te), (A.9) 



^ '2 '.'?>' 



idt2dt3dt4 



'2 

xQ(2'2)(t|4,fi,f2)AR(ii -f3)AR(i2 -f4)Q^^'^n*3j4|ii,4), 



(A.10) 
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t'2 



= - / dtidt2dt3dt4 



xQ(2'2)(ii,i2|ti,t2)AR(ii - t3)AR(i2 - h)Q^^'^\t3,U\t[,t'2), (A.ll) 

The general diagram rule is, match propagator outputs to vertex inputs, and vice 
versa, and integrate over all internal time arguments. By Mayer's first theorem 
[8l[31], series for response cumulants contain only connected diagrams. Coefficients at 
diagrams may be worked out following Vasil'ev [8]. 

Appendix B. The narrow- band case as a resonance approximation to the 
broad-band case 

In this appendix we show that the narrow-band case may indeed be recovered by 
making the RWA in formulae for the broad-band case. This demonstrates consistency 
of the RWA with response transformations. Although as expected, this result is not 
automatic. For simplicity, we consider the case when the RWA is applied to all modes. 
In other words, we look at the change from Af = to Af = N. This will allow us 
to reuse notation of section 13.21 with minimal amendments. Extension to the RWA 
applied only to a subset of modes is straightforward. 

For purposes of this section, the narrow- band field operators are given by equation 
PT|) with M ^ N, and the broad-band ones — by equation (jl2p with Af — >■ 0. Unlike 
in section 13. 2[ both operators are now defined with the same set of modes. The 
condition of resonance, 

|w„ - ujo\ < ujQ, (B.l) 

now matters, so that E{x,t) and E^x,t) are indeed slow amplitudes. Using equation 
(jB.ip we can write simple relations between the nonresonant and resonant operators. 
In the resonance approximation, E{x, t) and E^{x, t) are, up to factors, the frequency- 
positive and negative parts of A{x,t): 

A{x,t) = A^+\x,t) + A'--\x,t), 

E{x,t) w iLJoe'"°*A^+\x,t), (B.2) 
El{x,t) « -iujne-"^'>*A^-\x,t). 
By definition, we assume a similar structure for the current and dipole operators, 
J{x,t) = + 

« -iwoe-'^°'D{x,t), (B.3) 
j'^-^{x,t)^iiJoe"^°*DHx,t). 

D{x,t) and D^x^t) are also assumed to be varying slowly. Similar relations are 
postulated for the external sources; they are found from (jB.2p . (IB. 31) replacing 
A{x,t) Ac{x,t), E{x,t) ^ Ec{x,t), J{x,t) -J> Jc{x,t), and D{x,t) Dc{x,t). 
For the combination occuring in the T-exponents (|43p and (j44p we have, 

Aj + AJ + ijo = ED^ + E^D'' + ED* + E'^b + EID + E'i D^, (B.4) 

where we use condensed notation By itself, this relation is exact (subject to 

E{x,t) and D{x,t) being exactly frequency-positive, and E^{x,t) and D^{x,t) — 
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exactly frequency- negative). Cancellation of the counter-rotating terms is due to the 
formula, (with f{t) and g{t) being arbitrary functions) 

^+00 



dtf(+\t)g(+\t) 



dtf^-\t)g^-Ht)^0, 



(B.5) 



which is in turn due to the standard formula for functions and their Fourier-images, 

'•+°° duj 



dtf{t)g{t) 



27r 



(B. 



In the T-exponents (j43ll and (j44| . equation (jB.4p survives only as an approximate 
relation, subject to the RWA. What prevents it from being exact under the time 
ordering are the theta- functions in equations pT)) and p^ . Clearly that (|B.4[) survives 
as an approximate relation under the orderings is just another form of the usual 
argument, that fast counter-rotating terms must average out in equations of motion 
for slow amplitudes. The narrow-band S-matrices (|81|) and (I82p thus indeed emerge 
as a resonance approximation to the broad- band ones (|43p and (^H) , as expected. 

Consider now the "fate" of the combination ?7+A_|_ —ri-A- in equation (H^ under 
response transformation (I58p . Disregarding the time orderings we have. 



/ dxdt 

iujo 



3c , (-) 

ro(-) 



{+) 



A 



(-)l 



3c (+) 



(+) 



A 



(-)l 



3 — ILJot 



+'7^ '[x,t) 



E+{x,t) 



El{x,t) 



e-"^°*jt\x,t) 



E-{x,t) + e 



v'+>ix,t) 



El{x,t) 



(B.7) 



Because of the exponential factors condensed notation (|34l) becomes inapplicable, 
forcing us to write the last expression in full notation. Again, by itself equation 
(|B.7[) is exact; it turns approximate under the orderings. Comparing equation (jB.7| 
to ([83|) leads us into defining. 



e^'^o' j^+\x,t) _ 
fx+{x,t) ^ — ■ z , lJ'+(x,t) 



IJ,-ix,t) 



lOJQ 

-iuot ■{ — ) 
3e 



h 
{x,t) 



li-(x,t) 



e-'"°*r (_). , \x.t) 

— W '{x,t) + 



In turn, comparing this to the first line of 
the auxiliary variables. 



V^+\x,t) 



3i'^\x,t) 



(B., 



dc{x,t) ^ ■ , d^{x,t) 



we find the RWA correspondences for 

e-*"o*ji"^(a;,i) 



e'^o*ry(+)(a:,t) 
^i[X,t) = : ,^ (x,t) 



With these definitions. 



(B.9) 



fi*E+ - fiEl 



d: E: 



E^ 



+ 



djE 



El 



(B.IO) 
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The last formula here prepares ground for the response transformation of Wick's 
theorem. By similar arguments, 



ee{x,t) = : , e^{x,t) = -- 



e^'^°X(+\x,t) e-*-o*^(-)(a.^i) 

V[X,t) — ■ , V (x,t) — 



(B.ll) 



and 



= v*D+-vDl+ '^^ +^ ^ + zZ, (B.12) 

where the relation between p.±{x,t) and iy{x,t), ec{x,t) is by definition given by the 
second line in equation 

These formulae exhibit a number of important consistencies. The connection 
between de{x,t) and je{x,t) coincides with that between D{x,t) and J{x,t) and ipso 
facto with that between Dc{x,t) and Jc{x,t). The property that Jc{x,t) and jc{x,t) 
occur as a sum thus neatly turns into the property that De(x,t) and de{x,t) occur 
as a sum. Similar consitencies exist among pairs E(x, t), A(x, t), Ee{x, t), Ao{x, t) and 
ec{x,t),ac{x,t). Furthermore, equations (|B.9p force jc{x,t) to be real, and rj{x,t) to 
be purely imaginary. This agrees with the remarks in the last paragraph of section 
15.21 The same applies to ae{x,t) and v{x,t) in equations (IB. lip . 

To conclude our argument we look at the reordering forms. Using equations (P7| 
and [28)) it is easy to see that, 

ArXx, x\ t-t')^\ [e-*'^«(*-*')GR(x, x', t-t')+ e*'"''(*-*')G^(a;, x' ,t - t')] . (B.13) 
Making use of equations (jB.9p and dropping nonresonant terms we find, 



Unlike (IB.4p and (|B.7p , this relation is approximate because of the theta- function in 
Ar and Gr . 

Extending equation (jB.14p to quadratic forms of derivatives in equation (j70l) and 
(|73p takes a bit of mathematical effort. In the genuine narrow-band case of section [SI 
functional variables v{x,t) and ec{x,t) are unconstrained. Whereas their namesakes 
introduced by equation (jB.lip are, up to a factor, the frequency-positive parts of 
other variables. To show, for instance, that equation (|73|) turns under the RWA into 
(|102p . we need a bridging concept of constrained derivative by the frequency-positive 
and negative parts of a function. We define it by the formula, with f{x,t) being an 
arbitrary functional variable, 

(T) 

(B.15) 



6f(±){x,t) 
so that. 



Sf{x,t) 



sf{x,t) sfMix,t)^ sf(-){x,ty ^^-^^^ 

This definition follows the observation that separation of the frequency-positive and 
negative parts of a function. 
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f{x,t)^f<^+\x,t) + f^-\x,t), (B.17) 
is in fact an orthogonal decomposition. Indeed, for any pair of functions, 

dt/W(x,i)[g(-)(x',t)]*=0. (B.18) 



This formula may be verified using equation (jB.6p V Furthermore, the operations 
are projections. The frequency-positive part of a frequency-positive function is this 
function, the frequency-positive part of a frequency-negative function is zero, and 
so on. Definition (jB.15[) naturally extends the said orthogonal decomposition to the 
derivative d/df{x,t). In particular, for any functional F{-), 

6Fif)^ J dxdtj^^6f{x,t) 



dxdt 



Algebraic manipulation of derivatives (|B.15I) reduces to using their definition and 
observing that 5/5f^^\x,t) and 5/5f'^^\x,t) are frequency-negative and frequency- 
positive, respectively. It is then straightforward to show that 

so that equation (|102p is indeed a resonance approximation to (1751) . All formulae for 
the narrow-band case may thus be obtained making the resonance approximation in 
the formulae for the broad-band case. 

Note that, strictly speaking, we have demonstrated that the RHS's of all relations 
in section[5]are the RWA to the RHS's of the corresponding relations in sectionlH This 
is obviously consistent with the assumptions one normally makes about their LHS's. 
For instance, it is assumed that equations (jB.2p for the interaction-picture operators 
extend to the Heisenberg operators, 

A{x,t) ^ A^+\x,t) + A^-\x,t), 

£{x,t) « tLUoe''^"'A'^+\x,t), (B.21) 

£''{x,t) w -iu;oe~''^°'A^-Hx,t). 

The accuracy with which this formula holds is exactly the accuracy with which 
equation (jB.4p holds under the time orderings. The same applies to equation (|B.3p 
generalised to the Heisenberg operators. 
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